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Physiologically Based Pharmacokinetic Models 
And Biologically Based Pharmacodynamic 

Models 


LIBRARY OF PUBLISHED MODELS 


'Originally posted at PBPK-L Source Code Listing, ’ 

'URL: http://wvvw.navy.al.wpafb.af.niil/bew.htmi 1 
'and http://brandywine.members/plowchdr/ ’ 

'DISCLAIMER: By using these PBPK source codes you, the user,' 

' implicitly agrees with the conditions and rules stated in this ’ 

’ public domain source code disclaimer. If you have not read ' 

' this yet, please read it now. If you do not agree with the' 

' conditions and rules stated in it, please dispose of all contents r 
' of the PBPK source file and all related files now in your' 

' possession. The PBPK computer modeling programs are placed here by' 

' their authors for the convenience of others interested in PBPK ' 

' modeling. No representation of their suitability or correctness is’ 

’ implied. Please understand that PBPK models are complex entities * 

' whose value is minimal when taken out of context. With that' 

' understanding, you are welcome to use this code for teaching and ' 

‘ further development of the ideas that they represent. The' 

’ information and/or opinions presented here are those of the ' 

’ authors and are not endorsed by my previous or present employers, the US Air Force, US Navy, US ' 
' Army, TN & Associates, Inc., Department of Defense, or the US Government.' 

’ The US Air Force, US Navy, US Army, TN & A, Inc.,' 

’ Department of Defense, or the US Government ’ 

’ and the authors of this PBPK source code accept no responsibility ’ 

* for damages resulting from the use of the code and make no' 

' warranty or representation, either express or implied, including ’ 

' but not limited to, any implied warranty of merchantability or ’ 

’ fitness for a particular purpose. PBPK source code is provided *as’ 

’ is*, and you, its user, assume all risks when using it.' 


Go here to leant more about ACSL (simulation language) 



MODEL 1: LACTATIONAL TRANSFER OF TETRACHLOROETHYIJiNE AND CANCER 
RISK ASSESSMENT FOR INFANT 


DOWNLOAD MODEL FILES WRITTEN IN ACSL 


http: //members. spree. com/tnatox/PBPK/PBPK. htm 


11/15/00 

PM3006569703 


Source: https://www.industrydocuments.ucsf.edu/docs/fsdj0001 
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For LACTMAN.CSL file click here 
For LACTMAN.CMD file click here 

'Key words: Physiologically based phaimacokineric model, PBPK,' 
1 Tetrachloroethylene, PERC, PCE, Cancer risk assessment, ’ 

' Breast-fed infant, Toxicology, ACSL.' 


’ Description of this program has been published: Byczkowski,J.Z.,' 

' and FisherJ.W. (1995) A computer program linking physiologically ’ 
' based pharmacokinetic model with cancer risk assessment for ’ 

1 breast-fed infants. Computer Methods Progr. Biomed. 46, 155-163. ' 

1 Reprints of this article are indicated as AL/OE-TR-1995-0014 by ’ 
Armstrong Aerospace Medical Research Laboratory, Wright Patterson * 

1 Air Force Base, OH. ’ 

’ Results of simulation have been published: ByczkowskiJ.Z. (1996) 1 
1 A linked pharmacokinetic model and cancer risk assessment for ’ 

’ breast-fed infants. Drug Information J. 30,401-412. 1 
1 Reprints of this article are indicated as AL/OE-TR-1996-0007 by ' 

' Armstrong Aerospace Medical Research Laboratory, Wright Patterson ’ 

' Air Force Base, OH. ' 



MODEL 2: LINEAR MODEL, FOR CADMIUM INHALATION 


DOWNLOAD MODEL FILES WRITTEN IN ACSL 


For KNOB.CSL file click here 
'Key words: Linear pharmacokinetic model, PK, ’ 
' Cadmium, Cd, Inhalation, Risk assessment, 1 
’ Toxicology, ACSL. ’ 


’An ACSL version of the Kjellstrom andNordberg's PK linear model for’ 
'inhalation exposure to cadmium. ’ 

'Description of this program has been published: Byczkowski,J.Z.:' 
'Phamiacokinetic/Dynainic Modeling. In: Standard Handbook of' 
'Environmental Science, Health and Technology, Chapter 8.6 ' 

'(Lehr, J. Ed), McGraw-Hill, New York, 2000. ’ 



MODEL 3: LIPID PEROXIDATION IN LIVER SLICES 


DOWNLOAD MODEL FILES WRITTEN IN ACSL 


For LPOXL1V.CSL file click here 


http: //members, spree. com/tnatox/PBPK/PB PK. htm 


11/15/00 

PM3006569704 


Source: https://www.industrydocuments.ucsf.edu/docs/fsdj0001 
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For LPOXUV.CMD file click here 

'Key words: Biologically based pharmacodynamic model, BBPD, Liver' 
' in vitro, Lipid peroxidation. Toxicology, tert-Butyl hydroperoxide,' 

' Bromotrichloromethane, Mouse, Rat, ACSL.' 


!)’ Description of this program has been published: ByczkowskiJ.Z., 
’ Channel,S.R., Pravecek,T.L., and Miller,C.R. (1996) ' 
Mathematical model for chemically induced hpid peroxidation in ’ 
precision-cut liver slices: Computer simulation and experimental' 
calibration. Computer Methods Progr. Biomed. 50, 73-84. ’ 

Reprints of his article are indicated as AL/OE-TR-1996-0081 by ' 
Armstrong Aerospace Medical Research Laboratory, Wright Patterson ' 
Air Force Base, OH.' 

Results of simulation and experimental details have been' 
published: ByczkowskiJ.Z., Channel,S.R., and Pravecek.L. (1995).' 
Development and experimental calibration of the mathematical ’ 
model for lipid peroxidation in mouse liver slices. Reprints of' 
this article are indicated as Technical Report, AL/OE-TR-1995-0179' 
by Armstrong Aerospace Medical Research Laboratory, Wright' 
Patterson Air Force Base, OH. ’ 



MODEL 4: TCE - INITIATED LIPID PEROXIDATION IN LIVER SLICES 


DOWNLOAD MODEL FILE WRITTEN IN ACSL 


For LPOTCE.CSL file click here 

'Key words: Biologically based pharmacodynamic model, BBPD, Liver, ’ 
' in vitro, Lipid peroxidation. Toxicology, Trichloroethylene,’ 

' TCE, Mouse, TBARS, Thiobarbituric acid reactive substances, ACSL' 


’ Description of this program has been published: ByczkowskiJ.Z.,' 
' Channel,S.R., and Miller,C.R. (1999)' 

' A Biologically Based Pharmacodynamic Model for Lipid Peroxidation' 

' Stimulated by Trichloroethylene In Vitro.' 

' Journal of Biochemical and Molecular Toxicology. 13, 205-214.' 



EXAMPLE ACSL FILE OF BBPD MODEL 3 WRITTEN IN ACSL 


PROGRAM: LIPID PEROXIDATION 


'Version 2 *SQRT* (original model; Byczkowski et al. Comp. Meth' 

'Progr. Biomed. 50,73,1996), based on mathematical description’ 

'by Tappel et al. Free. Rad. Biol. Med. 7, 361, 1989. ’ 

'THE PURPOSE OF THIS PROGRAM IS TO SIMULATE PRODUC TION OF TBARS' 
•FROM PUFA IN MOUSE LIVER SLICES in vitro. **DEFAULT FOR TCE***’ 
'DESCRIPTION: Byczkowski et al. J. Biochem Mol. Toxicol. 13, 205, 1999.' 


http://members.spree.com/tnatox/PBPK/PBPK.htm 


11/15/00 

PM3006569705 


Source: https://www.industrydocuments.ucsf.edu/docs/fsdj0001 
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INITIAL 

’PEROXIDATION PARAMETERS FOR TCE IN LIVER [0.1 g of liver]' 
CONSTANT IND 1 = 1.0 S’inducer 1 [mcM TCE/0.1 g liv] ’ 

CONSTANT PUE - 7.0 $'Conc.PUFA[microMol/0.1 g liver]' 

CONSTANT PXZLUF= 12.0 S'peroxidizability of PUFA-L[l/h]' 

CONSTANT PXZHUF= 24.0 S'peroxidizability of PUFA-H[l/h]’ 

CONSTANT ANOX1 = 0.0037 S’Vit.E antiox[microMol/0.1g liv]' 

CONSTANT ACT1 = 0.0003 $'activatorl[mcM cytP450/0.1gliv]' 

CONSTANT GSH = 0.6 S'glutathione [mcroM/0.1 g liver]' 

CONSTANT PXTTBA= 0.1 S'yield ofTBARS/Mol hydroperoxide' 

CONSTANT LPUF = 3. $'LA-derivative PUFA[mcMo]/0.1 g]’ 

CONSTANT HPUF = 4. S'HA-derivative PUFA[mcMol/0.1 g]' 

CONSTANT EFANOl — 1. $’effectiveness of Vit.E ’ 

CONSTANT ANOXUF = 2. S'antioxidant use factor [/mcmol]’ 

CONSTANT ANOX2 = 0. $'non-Vit.E antiox[mcM/0.1 g liv]' 

CONSTANT ANOX3 = 0. $’added antioxidant[mcM/0.1g liv]' 

CONSTANT ACACT1 — 1. S'activity of activator 1 ' 

CONSTANT ACTDGF = 0.0014 S'activator degradation fctr /mem' 

CONSTANT PTIND1 = 250. S'potency of inducer 1 [1/mcmol] ' 

CONSTANT INDLF = 0.0001 S’inducer loss factor [1/h] ’ 

CONSTANT PXRATE — 0.00029 $'peroxidation rate [mcm/mcm] ’ 

CONSTANT AUTOXF = 0.00013 S'autooxidation factor [1/h]' 

CONSTANT GPENZA = 1. $'glutathione peroxidase[/mcM] 1 
CONSTANT PXREDF = 0.17 $'hydroperoxide reduction fctr /If 
CONSTANT TSTOP = 2. $'end of simulation [h] ' 

CONSTANT ACR = 0.025 $'Activator degradation rate[l/h]’ 

CONSTANT PHYSPX = 0. S'physiologieal levels of Hperox' 

CONSTANT ANRBG = 0. S'antioxidant regenerated in situ’ 

CONSTANT BCKGD = 0.0 S’TBARS in control [mcm/0.1 g] ’ 

!PARAME TERS FOR EXPOSURE TO SECOND CHEMICAL PROOXIDANT 
CONSTANT EFAN02 =- 0. S'effectiveness of non-Vit.E a-o ' 

CONSTANT EFAN03 = 0. S'effectiveness of added a-o ' 

CONSTANT ACT2 = 0. S'activator 2 [mcM/0.1 g liver]' 

CONSTANT ACACT2 = 0. S'activity of activator 2 ’ 

CONSTANT IND2 = 0. $'inducer2[mcM chemical/0. lg liv]' 

CONSTANT PTIND2 = 0. S’potency of inducer 2 [1/mcmol]' 

’TIMING COMMANDS' 

CONSTANT POINTS = 500. 

CINT=TSTOP/POINTS 

' INITIALIZATION' 

’RESETS INITIAL CONDITIONS BEFORE PEROXIDATION ’ 

CONSTANT ACTLOS=0. S'activator loss ' 

CONSTANT AUTOX =0. S'autooxidation ' 

CONSTANT PXREDG=0. S'Hperoxides red. by GSH peroxidase' 

CONSTANT PXLUF =0.00199 S'L-Hydroperoxides formed 1 
CONSTANT PXHUF -0.0053 $'H-Hydroperoxides formed ' 


http://members.spree.com/tnatox/PBPK7PBPK.htm 


11/15/00 

PM3006569706 


Source: https://www.industrydocuments.ucsf.edu/docs/fsdj0001 
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CONSTANT PXREM =0. S'accumulated remaining Hperoxides ’ 
CONSTANT TPX =0. $’accumlted total Hperoxides formed' 
CONSTANT ILR =0. S'rate of inducer loss = 0. ’ 

END $’End of Initial' 

DYNAMIC 

ALGORITHM IALG = 2 $'Gear stiff method ’ 

DERIVATIVE 

PROCEDURAL 

IF (LPUFRE. LEO) LPUFRE = 0. 

IF (HPUFRE.LE.0) 1IPUFRE = 0. 

IF (ACTLOS.GE. ACTEF) ACTLOS = ACTEF 
IF (ANOXRE.LE. l.e-10) ANOXRE = l.e-10 
IF (GSHREM.LE.O) GSHREM = 0. 

IF (INDLOS. G E. INDEF) INDLOS = TNDEF 
IF (PXREDA.GE.TPX) PXREDA = TPX 


END $'End of procedural ’ 


'#s Coixespond to equations in Tappel et al. (1989) and' 

’ Byczkowski et al. (1996)' 

!VARIABLES ARE CONVERTED TO CONCENTRATIONS [micromole/0.1 g liver] 


' 1. Remaining polyunsaturated fatty acids ’ 

LPUFRE = LPUF - PXLUFA - AUTOXA/2 
HPUFRE = HPUF - PXHUFA - AUTOXA/2 

'2. Effective activator' 

ACTEF = ACTDACACT1 + ACT2*ACACT2 
'3. Activator loss ’ 

ACTLOS = ACTEF* ACTDGF*TPX 
’4. Remaining activator' 

ACTREM = (ACTEF - ACTLOS)*exp(-ACR*INDREM*t) 
'5. Effective inducer’ 

INDEF = SQRT(IND 1 *PTTND 1) + IND2*PTIND2 


http://members.spree.com/tnatox/PBPK/PBPK.htm 


11/15/00 

PM3006569707 


Source: https://www.industrydocuments.ucsf.edu/docs/fsdj0001 
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’6. Remaining inducer' 

INDREM = INDEF - INDLOS 
'6.a. Inducer loss rate ' 

ILR = INDEF*INDLF 
'7. Activated inducer' 

ACTIND = INDREM * ACTREM 
'8. Effective antioxidant' 

ANOXEF = ANOX1 *EFAN01 + AN0X2*EFAN02 + AN0X3 ;i EFAN03 
'9. Remaining antioxidant 1 

ANOXRE = ANOXEF - ANOXEF*TPX*ANOXUF + ANREG*exp(- ANOXR*t) 

'lO.Hydroperoxides formed by action of activated inducer on' 

’ PUFA' 

PXLUF = LP UFRE *PXZLU F * ACTIND *PXR ATE/A NOXRE 
PXHUF = HPUFRE*PXZHUF*ACTIND*PXRATE/ANOXRE 

’ 12.Autooxidation ’ 

AUTOX = (LPUFRE + HPUFRE)*AUTOXF*TPX/ANOXRE 
1 14. Accumulated total hydroperoxides formed 1 
TPX = AUTOXA + PXLUF A + PXHUFA + PHYSPX 
'15.Remaining glutathione' 

GSHREM = GSH - PXREDA 

' 16.Hydroperoxides reduced by glutathione peroxidase ' 

PXREDG = PXREM*GPENZA*GSHREM*PXREDF 
'18. Accumulated remaining hydroperoxides ’ 

PXREM = TPX - PXREDA 

' 19.Amount of TB ARS from accumulated re mainin g hydroperoxides' 

TBARS = PXREM *PXTTB A + BCKGD 

http://members.spree.com/tnatox/PBPK/PBPK.htm 


11/15/00 

PM3006569708 


Source: https://www.industrydocuments.ucsf.edu/docs/fsdj0001 
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'6b. Inducer lost over time ’ 
INDLOS = INTEG(ILR,0.) 


' 11 .Accumulated hydroperoxides formed by action of activated' 
' inducer on PUFA ' 

PXLUFA - INTEG(PXLUF, 0.) 

PXHUFA = INTEG(PXHUF, 0.) 

'13.Accumulated autooxidation' 

AUTOXA = INTEG(AUTOX, 0.) 


'17. Accumulated hydroperoxides reduced by glutathione ' 

' peroxidase' 

PXREDA - INTEG(PXREDG, 0.) 

f _ r 

TERMT(T. GE.TSTOP) .^'Te rmin ation at TSTOP ’ 

END $'End of derivative 1 
END $’End of dynamic ’ 

'AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA' 

TERMINAL S’Resets parameters to initial values' 

ACTREM = ACTEF 
INDREM - INDEF 
ANOXRE = ANOXEF 
END $'End of terminal' 

f AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA' 

END $'End of program ’ 


Please send comments to the author E-mail address: 
j ANUSZ. B Y CZKOWS K1 ®\) S A. NET 

i i 


El 


http://members.spree.com/tnatox/PBPK/PBPK.htm 


11/15/00 

PM3006569709 


Source: https://www.industrydocuments.ucsf.edu/docs/fsdj0001 










